AC Josephson Effect Induced by Spin Injection 
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Pure spin currents can be injected and detected in conductors via ferromagnetic contacts. We 
consider the case when the conductors become superconducting. A DC pure spin current flowing 
in one superconducting wire towards another superconductor via a ferromagnet contact induces AC 
voltage oscillations caused by Josephson tunneling of condensate electrons. Quasiparticles simul- 
taneously counterflow resulting in zero total electric current through the contact. The Josephson 
oscillations can be accompanied by Carlson- Goldman collective modes leading to a resonance in the 
voltage oscillation amplitude. 

PACS numbers: 72.25.Dc, 71.70.Ej, 73.40.Lq 



I. INTRODUCTION 



Electric and spin transport near ferromagnct- 
paramagnct interfaces received a large attention boost 
with the discovery of the giant magnetoresistancc cffcctA 
and the subsequent developments in magnetoelectron- 
ics and spintronics. Aronov 2 -, and later Johnson and 
Silsbee^ theoretically predicted that an electric current 
through such an interface leads to an accumulation of 
nonequilibrium spin polarization with an accompanying 
spin current in the paramagnetic metal. A reverse ef- 
fect also takes place, a pure spin current from the nor- 
mal metal gives rise to an electric potential difference 
in the ferromagnet. The physics of these phenomena 
is quite simple. A sufficiently large difference of con- 
ductivities of spin-up and spin-down electrons in ferro- 
magnets induces spin polarization of the electric cur- 
rent therein. Spin polarized currents passing through 
ferromagnet-paramagnct boundaries result in an accu- 
mulation of nonequilibrium magnetization near the in- 
terface. Both the spin injection and detection of this 
spin polarization has been experimentally demonstrated 
in Refs. 0-0 in systems containing two or more junctions 
of thin normal metal wires with ferromagnets. One of 
them acts as a spin injector, while the other is a detec- 
tor, where the voltage created by diffusing spins can be 
measured. Related spin-polarized transport phenomena 
have been investigated in many spintronic applications, 
such as giant magnetoresistancei, spin Hall effects^, cur- 
rent induced magnetization dynamics^, spin-pumping^, 
and spin caloritronics^. 

In the case of superconducting systems, spin injec- 
tion and detection within a nonlocal setup similar to 
the one studied in Ref. 0-[l was investigated both 
theoretically-ii~— and experimentally^. These studies 
have been focused on DC transport. They revealed a 
strong renormalization of spin-related transport param- 
eters as compared to normal systems. These changes 
were mostly caused by the modified density of states in 
a superconductor. Beyond such quasi-particle transport 
properties, the macroscopic coherent state of the super- 
conducting condensate can give rise to a quite different 



transport phenomenon associated with the spin-polarized 
transport. 

Below we will consider an AC effect produced by a DC 
spin current towards a thin ferromagnetic contact. The 
DC potential induced by this polarization flux gives rise 
to an AC electric current of condensate electrons. Since 
in the considered experimental setup the total current 
of the superconducting and normal components must be 
zero, the AC condensate oscillations result in an AC po- 
tential difference between the opposite sides of the con- 
tact. A schematic of a possible experimental setup is 
shown in Fig. 1. A current is passed from a ferromagnet 
to a normal metal generating an associated spin accumu- 
lation and spin current therein. In the non-local geom- 
etry, this spin accumulation also diffuses transversely in 
a contacted normal metal towards another normal metal 
reservoir via a ferromagnet contact. The non-local poten- 
tial V increases with the injected DC current / and the 
nonlocal resistance R n i = V/ 1 describes the spin trans- 
port properties in the device. We will demonstrate that 
when the normal metals become superconducting, R n i 
acquires an AC component in addition to the DC com- 




FIG. 1: A possible setup for observation of Josephson volt- 
age oscillations. A spin current (dashed arrows) is injected 
by passing a DC current (solid arrow) from the ferromagnetic 
contact (brown). The DC spin current through the ferromag- 
netic contact between the right (R) and left (L) supercon- 
ducting electrodes induces periodic oscillations of their elec- 
tric potential difference V. 



2 



ponent in the normal state. 

The outline of this paper is as follows. A model sys- 
tem used in our calculation is described in Sec. III. Also 
in this section we present a simple calculation of the 
AC voltage oscillations assuming a local thermodynamic 
equilibrium between quasiparticles and the condensate. 
A microscopic analysis based on coupled kinetic equa- 
tions for the superconducting order parameter and the 
quasiparticlc distribution function will be given in Sec. 
III. A discussion of results will be presented in Sec. IV. 



II. AC VOLTAGE OSCILLATIONS IN A LOCAL 
THERMODYNAMIC EQUILIBRIUM 

Our model system consists of two superconducting 
wires in contact via a spin-active barrier. We consider 
this contact to be weak in the form of a thin ferromag- 
netic layer with, if necessary, additional insulating layers. 
Such a barrier can be characterized by two resistances 
Rf and R± corresponding to two spin eigenstatcs. We 
assume that a non-equilibrium spin polarization is cre- 
ated in the left wire (see Fig.l), either by spin injection, 
as shown in Fig. 1, or by other means. Moreover, we as- 
sume that the electron's energy relaxation is faster than 
their spin relaxation, so that up and down spin distri- 
butions can be characterized by the respective chemical 
potentials li^ and ^ , resulting in the spin accumulation 

potential 5fi s = (fj,^ — /x^J. In the right (R) wire Sfi s 

is much smaller, if the spin relaxation is faster than the 
influx of polarization from the left reservoir through the 
ferromagnetic contact. This is satisfied when the contact 
resistance is much larger than the resistance of the wire of 
the length l s = \/Dt s , where D is the diffusion constant 
and t s is the spin relaxation rate. This is true in many 
practical cases, in particular, in the systems studied in 
Ref. [BHll We, therefore, simplify our model assuming 
^ = fj.^ = fx — eV/2 and [t^ + fj,^ = 2ll + eV, where \x 
is the equilibrium chemical potential and V is the charge 
potential difference between two wires. 

With these definitions the electric current through the 
contact in the normal state is 



In _V , Sn, 



(i) 



where the inverse charge and spin resistances are given by 
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R, and i?„ 



t iii aim n s — iu. — ii, , respectively. 
In an open circuit, electro-neutrality requires I = 0, and 
Eq. (TTJ gives V = —Sfj, s R c /2eR s = —5ii s P/2e, where 
P = is the spin current polarization of the 

contact. This is just the voltage induced by the spin 
current trough the contact, as it has been experimentally 
demonstrated in Ref. 0-0- 

Let us now consider this situation for superconduct- 
ing wires. We assume that 28ll s <C ksT c , so that 
the nonequilibrium spin polarization does not cause 
depairingi^. The difference between Cooper pair energies 



on opposite sides of the contact is 2eV. This potential dif- 
ference gives rise to the AC Josephson current Ij of con- 
densed electrons. In addition, electro-neutrality causes 
an oppositely directed current /„ of quasi-particles, so 
that the total electric current is zero. This results in DC 
and AC voltage differences between the left and right 
superconductors that we will now compute. 

The simplest approach to this problem is based on the 
assumption that in the vicinity of the critical temper- 
ature T c — T -C T c , the quasiparticle current remains 
expressed by Eq. (fT]). We will discuss in the next sec- 
tion in which regime this approach is valid. Denoting 
the phase difference of the order parameters between the 
left and right wires as 0, and taking into account that 
d<fi/dt = —2eV/h, electro-neutrality I n + Ij = and Eq. 
(fT]) dictates 



I r sin < 



h d(p Sfi s 



2eR r dt 2eR s 



0. 



(2) 



where I c is the critical Josephson current. It is easy to 
see that when I c <C 5/i s /2eR s the Josephson current is 
dominated by harmonic oscillations with the frequency 
ui = PS^ s /h = 2eVo/h. Hence, the voltage induced by 
the spin current is 



V(t) = -V Q - I c R c smuit. 



(3) 



The DC component — Vq of this voltage is exactly the 
same as in the case of normal metals. Additionally, V(t) 
contains a term that oscillates with a frequency deter- 
mined by the DC (normal state) contribution of the non- 
local signal. The magnitude of the oscillating voltage can 
be estimated by noting that at temperatures close to T c 
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(4) 



where A is the superconducting gap and £ is a dimen- 
sionlcss coefficient that takes into account the depairing 
effect inside the ferromagnetic contact layer—, leading 
to exponential suppression of the Josephson current and, 
consequently to small values of £. It should be noted 
that according to Eqs. ([3]) and (j4|) the oscillation ampli- 
tude I C R C does not explicitly depend on the transmission 
coefficient, apart from the weak dependence through the 
depairing factor Q. 



III. NONEQUILIBRIUM EFFECTS AND 
COLLECTIVE MODES 



The above analysis was based on the assumption that 
near the critical temperature, the current carried by 
quasiparticles can be represented by the expression in 
the normal state Eq. ([1]), ignoring small corrections as- 
sociated with the gap in the quasiparticle spectrum. The 
small gap alone, however, does not justify this assump- 
tion. In particular, when quasiparticles are transmitted 
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between the left and right wires they may not be in the lo- 
cal thermal equilibrium with the respective condensates, 
that have been assumed when deriving Eq.(|2]). To take 
into account nonequilibrium effects, one needs to consider 
time dependent transport and relaxation of the quasi- 
particles. There are two physical effects that determine 
kinetics of quasiparticles in the superconducting wires. 
The first one is the so called^— charge (or branch) im- 
balance of electron and hole excitations. It is produced 
by quasiparticle tunneling between superconducting elec- 
trodes, leading to a quasiparticle distribution with a local 
chemical potential different from that of the condensate. 
This difference relaxes during a time much longer than 
the electron-phonon scattering time. Another effect is re- 
lated to condensate space-time oscillations. It dominates 
over the charge imbalance relaxation when u> is large 
enough. We will demonstrate that the spin injection then 
enables detection of collective condcnsatc-quasiparticle 
modes, Carlson-Goldman modest which are character- 
ized by oppositely directed oscillations of condensate and 
normal fluids. There is an important difference with re- 
spect to the usual Josephson effect, since our device re- 
quires no net current 1 = 0. The usual Josephson effect 
does not couple to Carlson-Goldman modes and is not 
reduced at low temperatures, T — > 0. In contrast, to 
provide a counterflow we need excitations that vanish at 
low temperatures. The coupling to the collective modes 
is enabled by a spin-driven battery effect induced by the 
spin injection. 

Let us now detail the calculations. Assuming a small 
deviation from equilibrium we employ the linearized time 
dependent kinetic and Ginzburg-Landau equations in the 
diffusive regime^ 8 -, when the elastic mean free path is 
much less than the superconductor's coherence length, 
as well as other relevant length scales. In this case, 
the isotropic quasiparticle distribution function f a (E, t) , 
where a is the spin projection, depends only on the en- 
ergy E and time t. Within the linear theory the sin- 
glet condensate couples to the spin-independent part 
f(E, t) = (ff(E, t) + f±(E, t))/2 of the distribution func- 
tion. Therefore, after ignoring small terms (S/j, s /kBT) 2 , 
the unperturbed spin-independent distributions takes the 
form of Fermi equilibrium functions f (E, t) of the left 
and right wires with respective electrochemical poten- 
tials eV/2 + \x and — eV/2 + \i. In its turn, the corre- 
sponding gap functions of unperturbed condensates are 
Aexp(i(/)/2 — 2i/j,t) and Aexp(— itfi/2 — 2i[it). It is easy 
to see that in this unperturbed state the spin indepen- 
dent contribution to the quasiparticle current through 
the contact is given by the first term in the right-hand 
side of Eq. (jTJ . Taking into account above condensate 
functions one can easy obtain Eqs. @ and ©. In the 
perturbed state we have f(E, r, t) = fo(E, t)+6f(E, r, t) 
(we will skip here and below the labels L and R). Since 
the perturbation violates the electron-hole symmetry, it 
gives rise to a spatially dependent potential <^(r, t) near 
the contact. Also, a correction to the order parameter 
<5A(r, t) appears. In order to simplify the further analy- 



sis, we assume that tu <^ A and 1/te <C A, where 1/te 
is the electron-phonon relaxation rate. Besides that, the 
critical supercurrent I c is taken small enough, so that the 
time dependence of all functions is dominated by har- 
monic oscillations. Accordingly, we introduce the time 
Fourier components 5/^(r,i?), <5A w (r) and ^(r). From 
Refs. [I8II2TI it follows that obeys the kinetic equation 

(-itoN! - DV 2 + 2AA 2 ) 5f u + 

iuiNtfo&fu - LoN 2 f SA UJ = I st , (5) 

where f = 1 /4k B T cosh 2 {E/2k B T) and I st is the 
electron-phonon scattering integral, whose explicit form 
can be found in Ref. GJ|2ll Furthermore, D = D(N? + 
N$), where Ni and N 2 are the spectral functions: 

2^1/2 = G 1/2 (E + w/2) T G 1/2 (E - w/2), (6) 

where d{E) = E/y/E 2 - A 2 and G 2 (E) = 
«A/V E 2 — A 2 , with to = lo + i/te- In its turn, the lin- 
earized Ginzburg-Landau equation takes the form 

- iu5A u + 8 ' fc f A TA / dEN 2 5f w = DV 2 5A U . (7) 
tt|A| J 

We will employ the above equations for the analysis of 
our model in two limiting cases of weak and strong en- 
ergy relaxation versus the Josephson frequency, t^w< 1 
and teu 1, corresponding to very different physi- 
cal situations. In the former case slow time variations 
of 5f may be ignored, so that the quasiparticle kinet- 
ics is dominated by the charge imbalance of electron and 
hole excitations. The deviation from the thermodynamic 
equilibrium decreases with increasing distance from the 
contact on the characteristic length scale ^/Dtr, where 
tr = AksTcTE/^A is the charge imbalance relaxation 
time, that is much longer than te- In the opposite high- 
frequency regime inelastic collision processes are not im- 
portant, because the quasiparticle distribution oscillates 
fast. Therefore, one can neglect 1/te and I s t in Eqs. 
(|5I7[) . In this case, since Josephson oscillations of the 
condensate take place at zero total current, they strongly 
couple to Carlson-Goldman modes. Therefore, one can 
expect such modes to be excited near the contact and 
propagate along the left and right wires. 

In both low-frequency and high-frequency regimes, us- 
ing Eqs . (|5l7p with a reduced form of I s t from Refs. [2l] 
and [l8| and taking into account the zero electric current 
condition, one arrives to the equation for the potential 

k 2 (w)^=V 2 ^, (8) 
where k 2 (uj) = 1/Dtr at teoj <C 1 and 

C 2 S K 2 (LU) = -LO 2 - iTTL0A 2 /4k B T (9) 

at teoj S> 1, where the sound velocity c s = \/2DA. Eq. 
(j8]) is well known. At te^o <C 1 it describes the charge 
imbalance relaxatio n 17 ' 18 , while in the opposite limit it 
gives the dispersion of Carlson- Goldman modest. 
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For our geometry, when k -1 is much larger than the 
width and thickness of the wire, <p u depends only on the 
coordinate x along the wire. Then, at teuj < 1, ^ 
exponentially decreases with increasing distance from the 
contact, while at tew ^> 1 it shows decaying oscillations. 
We assume that the left and right wires are of the same 
length L. Since the system is symmetric with respect 
to x —> —x, the oscillating part of the electrochemical 
potential is — V u /2 + (p u (x) at x > and V^/2 — <^ w (— x) 
at :r < 0, where denotes the Fourier component of 
V(t). The solution of Eq. © has the form 



-k{uj)x 



with the boundary conditions V a; ^ Cd (zbL) = and 
V u -2<p u (0) 



= AaV x ip u (0), 



(10) 



(11) 



where A is the wire cross-section area and a is the normal 
state conductivity. These boundary conditions provide a 
zero electric current of quasiparticles at the wire ends and 
the current equal to the injected one at x = 0. From Eqs. 
([8|). (fT0)) - ([TT]) one obtains a periodic part of the injected 
current 



^-2^(0) 
Rc 



v., 



R c + 2R U 



Reff(u) ' 



(12) 



where 



Rut = 



1 1 



-2kL 



Aan(uj) 1 — e 



-IkL 



,Re(/t)>0. (13) 



Hence, the result is a renormalization of R c in Eq. ([TJ, 
such that R c — s- R c + 2R W . To find the voltage V^, 
the quasiparticlc current (|12[) must be equated with the 
Josephson current. By this way we obtain a new expres- 
sion for a time dependent part of V, instead of the second 
term in the right-hand side of Eq. ([3]): 



V+V 



-7 c [cos LutlmR e ff(uj) + sin u)tRcR e f j (u>)} (14) 



IV. DISCUSSION 

Let us analyse above results in some limiting cases. 
Since k(uj) — > oo in both cases of high frequencies u> — > oo 
and strong energy relaxation te — > 0, it follows from Eq. 
([TBI that Rw ->• 0. We thus obtain Eq. flBJ, that is an 
expected result, because in these limits a deviation from 
equilibrium is small. On the other hand, the noncqui- 
librium effect of quasiparticle's kinetics becomes strong 



when R c < 2R W . The assumed linearization condition, 
however, restricts this inequality. This condition can be 
expressed in the form C,\R W /R C \ <C 1. Therefore, the lin- 
ear theory allows R c < R w only at small (. It should be 
noted that, according to Eq. ([TB"]) , R w can be enhanced 
due to resonances of Josephson oscillations with collec- 
tive modes at IitikL = %n, if they are not overdamped 
(if RckL <C 1). R w also increases at small enough L, 
when 2 1 ft | L <C 1. In practice R w may be varied in quite 
wide range. In Al wires from Ref. [g Vq = 10 -6 V, re- 
sulting in lj- 1 ~ 0.3 • 10~ 9 s. Since w" 1 ~ t e - 10" 9 s^, 
a regime intermediate between charge imbalance relax- 
ation and generation of Carlson-Goldman modes will be 
realized, with kT 1 about several /im. Therefore, strong 
resonances in R w are not expected. One can evaluate 
R w ~ 50fi, that is much less than R c = 60051. Hence, 
in the considered parameter range Eq. (jB)) remains valid. 
In samples with higher polarizations P and at larger spin 
current through the contact the Josephson oscillation fre- 
quency is expected to be large enough to produce notice- 
able collective resonances of R e ff in Eq. (UJ|. 

The above calculations of Josephson voltage oscilla- 
tions have been restricted to A <C T c . At the larger gap 
the oscillation amplitude is expected to decrease, because 
less excitations are available to compensate the super- 
current through the contact. On the other hand, in this 
range one should take into account that besides quasipar- 
ticles the spin transport through the contact can be asso- 
ciated with triplet components of Cooper pair states that 
appear due to spin dependent tunneling and nonequilib- 
rium spin polarization of superconducting wires. Further 
studies are needed to understand the effect of such trans- 
port. 

In conclusion, we considered an AC Josephson effect 
induced by a DC spin current through the contact whose 
transmittance depends on the spin orientation of tunnel- 
ing electrons. The oscillations of the voltage across the 
contact at zero electric current have, in certain parameter 
range, harmonic time dependence with the frequency pro- 
portional to the spin current. The amplitude and phase 
of these oscillations depend on coupled kinetics of quasi- 
particles and condensate in superconducting wires. The 
corresponding calculations have been performed within 
linearized kinetic equations at temperature close to T c . 
We predict that at the high enough frequency the mea- 
sured AC voltage will show up the resonance structure 
associated with excitation of Carlson-Goldman modes. 
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